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ABSTRACT 

We use the Geneva Syclist isochrone models that include the effects of stellar 
rotation to investigate the role that rotation has on the resulting colour-magnitude 
diagram (CMD) of young and intermediate age clusters. We find that if a distribution 
of rotation velocities exists within the clusters, rotating stars will remain on the main 
sequence (MS) for longer, appearing to be younger than non-rotating stars within 
the same cluster. This results in an extended main sequence turn-off (eMSTO) that 
appears at young ages (~ 30 Myr) and lasts beyond 1 Gyr. If this eMSTO is interpreted 
as an age spread, the resulting age spread is proportional to the age of the cluster, i.e. 
young clusters (< 100 Myr) appear to have small age spreads (10s of Myr) whereas 
older clusters (~ 1 Gyr) appear to have much large spreads, up to a few hundred 
Myr. We compare the predicted spreads for a sample of rotation rates to observations 
of young and intermediate age clusters, and find a strong correlation between the 
measured ’age spread’ and the age of the cluster, in good agreement with models of 
stellar rotation. This suggests that the ’age spreads’ reported in the literature may 
simply be the result of a distribution of stellar rotation velocities within clusters. 
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1 INTRODUCTION 


The phenomenon of extended main sequence turn-offs 
(eMSTOs), where the MSTO in colour-magnitude dia¬ 
grams is broader than predicted by a single isochrone (in¬ 
cluding the effects of binarity and photometric errors), 
is a common feature of intermed iate age (1 — 2 Gyr; 
e.g. iMackev fc Brobv Nie lsen 20 o 3) and apparently even 
youn ger (~ 300 Myr; Milone et al.l l2015llCorrenti et ail 
l2015h massive stellar clusters. The origin o f thi s f eature 
is still highly debated in the literature (e.g. iLi et al. 1 1 2014 
IGoudfrooii et all 1 201,51 ). 

One interpretation, which disagrees with the long stand¬ 
ing paradigm of star clusters being simple stellar popula¬ 
tions (SSP), is that there is an age sprea d inside the c lus¬ 
ters of the o r der of ~ 200 — 700 Myr (e.g. G!oudfrooij_e^il 
20091.1201 lal fbl. I20l4 120151: iRubele et al.ll2013l:ICorrenti et al 


2014 ). As a scenario to form such clusters. IGoudfrooii et al 
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(l2011bl. 1201411 suggested that a massive cluster initially 
emerges as an SSP from a near-instantaneous star form¬ 
ing burst, which is then followed by a Gaussian-like ex¬ 
tended star formation episode. The second generation stars 
are formed out of a combination of the ejecta of evolved first 
generation stars and large amounts of accreted gas from the 
surroundings. The authors then invoke strong cluster dis¬ 
solution, due to tidal forces and cluster expansion, which 
preferentially removes 1st generation stars. By an age of 
1 — 2 Gyr, the clusters have removed the majority of their 
first generation stars, leaving nearly only the 2nd generation 
with an extended star-formation history (SFH). While po¬ 
tentially plausible for the intermediate age clusters, this sce¬ 
nario does not appear to be consistent with the inferred SFH 
of the younger cluster NGC 1856 (~ 300 Mvr: I Milone et all 
l2015h . where the r ole of cluster dissolutio n is expected to 
be very small ('e.g-. lBaumgardt et al.lfioi4 and two distinct 
peaks in the SFH are not seen. 

Several studies have aimed to search for signs for 
prolonged star formation histories or ongoing star forma- 
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tion in young massive clusters (e.g. BastiaB_gt_aL | 20131: 
ICabrera-Ziri et all 120141 : iNiederhofer et alT 201 Sal . Cab re r a- 
Ziri et al. in prep.), however no such extended or ongoing 
bursts have yet been found. Theses non-detections call into 
question the interpretation of age spreads as the cause of 
extended MSTOs. 

Apa rt from age spreads, other scenarios have been put 
forward. iBastian fe de Mink) ll2009ll were the first to sug¬ 
gest that the eMSTO feature may be due to stellar rotation. 
In their paper they computed the effect that rotation has 
on the surface colour and luminosity of stars that rotate 
at various velocities and were able to reproduce a broaden- 
ing in the turn-off s imilar to the observed ones. However, 
IBastian fe de Minkl l2009ll did not include the effect that 
rotation has on the interior of t he stars, e.g. a longe r MS 
life time due to internal mixing. IGirardi et al] (1201 ill crit¬ 
icised the proposed rotation scenario and showed that the 
prolonged lifetime of moderately rotating stars can compen¬ 
sate for the shift in colo ur and magnitu de, leading to a tight 
MSTO. In contrast. I Yang et all d2013l ') found in their study 
that the extended MSTO feature can, at least partially, 
be explained by stellar rotation, based on calculations per¬ 
formed with the Yale Rotational Stellar Evolutiona ry Code 
(YREC, iPinsonneault et al.lll989l : lYang fe Bill2007t) . 

More recently, studies have investigated the effects of 
rotation on the MSTO (and colour-magnitude diagrams 
mo re generally) based o n the Geneva rotat i ng stel lar mod¬ 
els llGeorgy et a.1 .1 Eli) 1 .3 al l . [Brand I fe Huand (l2015al l studied 
the CMD of the Praesepe and Hyades, two nearby Galactic 
open clusters (~800 Myr) that both have turn-offs which 
are not in agreement with a coeval stellar population. They 
showed that the spread in their turn-offs can be explained 
with rotating stars. The authors have exp anded thei r mod¬ 
els to older isochrones, up to 2.5 Gyrs ll Brandt fe Hu ansi 
l2015bl l. The models predict an extended MSTO whose width 
is a function of the cluster age with a maximum spread at 
around 1 — 1.5 Gyr. A comparison of their models with four 
extended MSTO clusters show that stars with the same age 
but different rotation rates and inclinations can reproduce 
the spread i n the observed turn-off. 

Finally, INiederhofer et ahl (l2015bl) have found a strong 
correlation between the width of the ’inferred age spread’ 
and the age of the cluster, with younger clusters showing 
smaller age spreads than older clusters below ~ 1.5 Gyr and 
the trend inverting beyond this age. The correlation coef¬ 
ficient between the age and FWHM is much stronger than 
between the FWHM and present day m ass or escape ve¬ 
l ocity (or even the ’corrected’ values from iGoudfrooii et al.1 
(120141) 1 that have been used to argue for a causal connection 
between the cluster properties and age spreads. It is inter¬ 
esting to note that the younger cluster NGC 1856 contin¬ 
ues this trend to lower ages, having a much smaller inferred 
spread (At = 100 — 130 Myr) than its older intermediate 
age clusters (At = 300 — 700 Myr). 

Stars are known to rotate since Galileo Galilei and Jo¬ 
hannes Fabricius observed the Sun’s spots around 1610. Very 
early after the understandin g of the stellar constitution by 
[Eddington! l)l916j . fl917l . ll918l l. the need to adapt the internal 
structure equations to the case of stellar rotation was studied 
(lMilnelll923l : Ivon Zeipellll924l : lEddingtonlll925l l. The estab¬ 
lishment of a large scale circulation in the interior o f rotating 
stars has been presented first by lEddingtonl dl929l l. and for- 



Figure 1. Four isochrones from the Syclist models. The 
isochrone of a rotating star (ui = 0.5) of age log t=8.5 closely 
resembles that of a non-rotating younger isochrone of an age log 
t=8.35 at the main-sequence turn off. Hence, a distribution of ro¬ 
tation rates within a young/intermediate age cluster has a similar 
effect as an age spread on the MSTO portion of the CMD. The 
points mark specific points in the evolution of the star and are 
described in the text. 


malised later by ISwee tl (Il950l) . This circulation (Eddington- 
Sweet or meridional circulation ), is responsible for a mix¬ 
ing inside the star that adds to the mixing by convective 
movements. The effects of rotation hav e been included in 
stellar evolution cod es since the 1970’s l|Kippenhahnl Il97(il : 
lEndal jsT'S ofia 1976). After the observations of enric hments 
in N I Lvubimkovlll984h and He dHerrero et al .1 ITooij) in OB 
stars, with larger enrichments in more massive stars, rota¬ 
tion was rapidly identified as the cause and started to be 
seen as a necessary ingredient in stellar models. While the 
numerical details of its implementation in the various stel¬ 
lar evolution codes differ in many ways, the general effects 
are always the same: 1) an internal mixing that modifies the 
chemical composition at the surface and increases the MS 
lifetime, 2) hydrostatic modifications on the structure. Both 
effects have strong consequences on the stellar evolutionary 
paths. 


2 SYCLIST MODELS 

In this st udy, we used the stellar population synthesis code 
Syclist llGeorgy et al.ll2014n . designed to be used with re¬ 
cent grids of models computed with the Geneva stellar e vo- 
lution code (lEkstrSm et all 120121 : iGeorgy et ah 2013a.2). 
Fo r the present work, the grids presented in Georgy et al.1 
l|2013alf are used and we adopt a metallicity of Z=0.006. 
They cover a mass range between 1.7 and 15 Mq and ini¬ 
tial rotation rates between 0 and 95% of the critical on<0. 
They are thus well designed for young and intermediate 


1 We use here the notation oj = S!/Q cr j, for the ratio of the actual 
surface angular velocity to the critical one. The critical angular 
velocity accounts for the effect of the deformation of the shape of 
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age clusters between ~ 20 Myr and ~ 1 Gyr, and allow 
one to consider an initial velocity d istribution. Bef o re the 
fina l implementat i on of t he grids in lEkstrom et, al.1 (|2012h 
and lGeorgy et ail (l2013al 'l was done, several tests with vari¬ 
ous prescriptions for rotation, including angular momentum 
transport and rotational mixing, were performed. The pre¬ 
scription which reproduces best observational features was 
finally selected. 

Syclist is built to compute interpolated models, 
isochrones, and synthetic clusters. By summing up clusters 
of various ages (with different weights), it is also able to 
compute stellar populations with different SFHs. The ex- 
act way stel lar populations are computed is described in 
iMevnet et al.l (120151 ). In Syclist, the colours and magni¬ 
tudes are obtained from the luminosity, surface gravity, ef- 
fective tempera t ure an d metallicity following the work of 
IWorthev fc Lee! (|201ll 'l. In this paper, we us e the obser- 
vation ally determined velocity distribution of iHuane et al.1 

d20inh . 

The effect of rotation on CMDs of clusters was al¬ 
ready studied in previous works. Due to the different ap¬ 
proaches followed so far, we summarise here briefly the 
m ain differences betwe e n the se studies and the present one. 
In lBastian fc de Minkl l|2009f ). the effects of rotation where 
parametrised from a test case, and then applied to non¬ 
rotating isochrones to mimic the effect of rotation. The 
main problem of this approach is that the modification of 
the lifetime of rotating stars is not accounted for (rotating 
stars have a longer lifetime than n on-rotating ones, at least 
in the framework o f the grids by iGeorev et al.l l2013al ) . In 
iGirardi et al.1 ([201ll j. rotation was accounted for the same 
way than we do here. However, differences can be noted 
about the choice of the diffusion coefficients used, as well as 
the convective overshooting. Moreov er, they focu sed on clus¬ 
ters with a higher age than here. In lYang et al.i d2013T ). the 
transport of angular momentum inside the star is treated as 
a purely diffusive process, contrarily to our approach where 
the advective term in the equation of angular momentum 
transport fseelZa h7il ll992i) is accounted for. 

2.1 Determining the Expected ’Age Spread’ 

As shown in iGeorgy et al. ! ll20l4 and references therein) 
rotation at young to intermediate ages causes heightened 
mixing within stars, allowing more hydrogen to enter the 
core, i.e. rotating stars stay o n the M S for longer than non¬ 
rotating stars (see also lGirardi et ahll201ll ). This manifests 
as an eMSTO that can mimic an age spread within clusters. 

In order to estimate what effect a rotation-induced 
MSTO spread would be if it was inferred to be an age spread, 
we performed the following tests, shown in Fig. [T] We first 
adopt an isochrone of a given age, e.g., log t = 8.5, and 
cj = 0.0 (i.e., a non-rotating isochrone), and compare it to 
an isochrone of the same age but for rotating stars with 
u> = 0.5 (which are labelled isochrones 1 and 2 in Fig. [T] 

the stellar surface produced by rotation (cf. jMaede^j^OOEj): 

8 GM 

27 Rl crit 

where R p cr j, is the polar radius at the critical velocity. 



Figure 2. The predicted and observed relation between the in¬ 
ferred age spread in young and intermediate age clusters and their 
age. The dotted (blue), the dashed (red) and the solid (black) lines 
show the prediction of the Syclist models if rotation-induced 
main sequence turn-off spreads are inferred to be age spreads 
based on the differences in magnitudes at three different points 
along the isochrone (see text). A key prediction of the rotational 
models is that the inferred age spread should increase with clus¬ 
ter age (until magnetic breaking becomes important). The filled 
points show observations of young and intermediate age clusters 
where age spreads have been inferred based on the MSTOs. The 
data are well reproduced by the predictions of rotation. 


respectively). We then over-plotted non-rotating isochrones 
for younger ages (labelled as isochrones 3 and 4 in Fig. |T} 
until the youngest isochrone that best reproduced the log 
t=8.5, uj = 0.5 isochrone was found. As shown in Fig. [X] in 
the MSTO region of the CMD, rotation manifests itself in a 
way that is nearly identical to that of an age spread. 

We carried out the above test in an automated way for a 
range of ages, from 50 Myr to 1 Gyr, and derived the inferred 
spread if the rotation induced broadening was interpreted as 
an age spread. For this we selected a specific point along the 
non-rotating and rotating isochrones and searched for the 
closest match in magnitude. In order to be more objective, 
we performed the analysis three times using three different 
points along the isochrones. The locations of these points are 
indicated as circles along isochrone 2 in Fig.Q] The first one 
is the brightest point (minimum My) after the stars have 
left the MS (blue dotted circle). As the second position in 
the CMD we selected the evolutionary MSTO (red dashed 
circle). As the last position we chose the point where the 
isochrone passes right above the local brightness maximum 
of the hook-like feature just before the MSTO (solid black 
circle). 

The results are shown as a dotted (blue), dashed (red) 
and solid (black) line in Fig. [5] All three choices of the spe¬ 
cific points along the isochrones yield similar results. The 
inferred age spread is proportional to the age of the clus¬ 
ter, with younger clusters showing smaller age spreads than 
older clusters. For clusters ^ 100 Myr, age spreads of a few 
10s of Myr would be expected, whereas for intermediate age 
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Figure 3. Upper panels: Hess diagrams centred on the MSTO 
features of two artificial clusters. Left panel: A cluster constructed 
from stars with different ages. The ages are from a Gaussian dis¬ 
tribution with a peak age of 250 Myr and a standard deviation 
of 50 Myr and are sampled at 11 equally spaced times. Super¬ 
imposed as black lines are two isochrones at 155 and 345 Myr. 
Right panel: A coeval cluster but with stars rotating at various 
velocities. The age of this popula tion is 300 M yr. The rotation ve¬ 
locity distribution is taken from iHuang et al.l j20ldl . The black 
lines are isochrones for non-rotating stars (lo = 0.0) and stars 
that rotate with half of their break up velocity (u = 0.5). Both 
diagrams are blurred with a Gaussian with a dispersion of 0.03 in 
colour and 0.02 in magnitude to account for photometric noise. 
Note that the colour scale in both panels is in logarithmic units. 
Lower panels: The same as the upper panels but now centred 
on the post-MS tracks of the cluster CMD. 


clusters (~ 1 Gyr), age spreads in excess of 300 Myr are 
expected. Note the similar relation of the observations. 

For older ages, rotation-induced spreads are expected 
to decrease at some point. The stars at the turn-off have 
sufficiently low masses to posses convective envelopes gen¬ 
erating magnetized winds. Through these winds the stars 
loose angular momentum which slows down the stars. This 
happens at a turn-off mass of ~1.3 Mq and is referred to as 
the break in the "Kraft curve" {Kraft 196^. _ 

Based on the results of Niederhofer et al] (l2015bh 
this appears to happen around ~ 1.5 Gyr, although 

iBrandt fe Huana (I2015br ) suggest that this will happen 
somewhat earlier, at an age of ~ 1 — 1.5 Gyr. 

We would like to mention here that we used for our anal¬ 
ysis rotating isochrones with an initial rotation rate u>=0.5 
which is not an extreme value compared to what is measured 
for B stars (e.g. IHuang et alJuOld f . The results that we ob¬ 
tained are valid as long as the initial rotation distribution 
roughly covers the interval from uj=0 to cj= 0.5. If the initial 
rotation velocities for specific reasons is peaked towards a 
single value, the eMSTO features cannot be explained by 
stellar rotation. 

Besides stellar rotation, also convective overshooting 
modifies the isochrones and could affect our results. How¬ 
ever, changing the overshooting in the models would only 
shift both isochrones, rotating and non-rotating, towards 


higher luminosities in the CMD, and therefore the effect 
on our results is minor. As the masses and compositions of 
the stars at the MSTO are nearly identical, variations of 
overshooting along the MSTO stars is not expected. 


2.2 Effect on the Post-Main Sequence Features 

An age spread within a cluster and a range of stellar rotation 
velocities will not only effect the MSTO but also the post- 
MS features in the CMD. To test how these two scenarios 
are reflected in the structure of a CMD, we created two syn¬ 
thetic clusters using the Syclist models and compared their 
relative stellar densities in colour-magnitude space(Hess di¬ 
agram). The first cluster is constructed out of non-rotating 
stars with a Gaussian-like age distribution that has a peak 
age at 250 Myr and a standard deviation of 50 Myr. The 
second cluster consists of a coeval stellar population with 
an age of 300 Myr, however, with stars rotating at various 
velocities. The rotation distribution wa s taken from the o b- 
servationally determined distribution of lHuang et al.1 lj201(t l . 
The upper panels in Fig. [3]show the B — V vs My Hess dia¬ 
gram of the two artificial clusters centred on the MSTO. The 
left panel shows the age spread cluster, whereas the cluster 
with the rotating stars is displayed in the right panel. The 
structure is slightly different in both models, however, the 
location of the turn-off itself (denoted by the top of the high 
density region in red) is nearly at the same magnitude and 
colour in both cases. Hence, age spreads and a rotation dis¬ 
tribution can lead to very similar MSTO features. 

The post-MS morphologies, however, are very different 
in the two scenarios, as shown in the bottom panels of Fig. [3] 
In the age spread cluster, the post-MS tracks extend over a 
large magnitude range (Fig. [3]left) whereas the evolved stars 
in the cluster with coeval rotating stars show a more com¬ 
pact morphology in colour-magnitude space (Fig. fright). 

The fact that the morphology of the post-MS dif¬ 
fers in the age spread and rotation model can help to 
discri mi nate between the two scenarios. We note that 
iBastian fc Silva-Villal (|2013h focussed their study largely on 
the post-MS features of NGC 1866 and 1856, and found rel¬ 
atively small age sp reads, i.e. o < 3 5 My r for both clusters. 
On th e other hand. ICorrenti et al.l d2015l f and lMilone et all 
d2015ll focussed on the MSTO of NGC 1856 and found evi¬ 
dence for a much more extended SFH, 90 — 130 Myr. This 
apparent contradiction is predicted if rotation is the under¬ 
lying cause of the observed extended MSTOs. 


3 COMPARISON WITH MEASUREMENTS 
FROM THE LITERATURE 

A number of papers have reported age spreads based 
on e MSTOs in y o ung a nd in t ermediate age c l usters 
(e.g., Miloneetahl l2009t 120151 : iGoudfrooii et alJ 12014 
ICorrenti et al.l 2015l h In Fig. [5] we over-plot the data from 
the literature on the expectations from the analysis dis¬ 
cussed in §EU The fill ed circles are in termediate age clus¬ 
ters from the sample of iGoudfrooii et al.l (I2014T ). where we 
have adopted their age and FWHM of the age distribu¬ 
tion for each cluster in the LMC (blue) and SMC (green). 
The filled (re d) square repr e sents NGC 1856 from the re¬ 
cent study of iMilone et al.1 (|2015l f and the (yellow) filled 
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squa re represents t he sam e cluster but the age spread is 
from ICorrenti et al.l I^ Ols HB- The filled (red) triangle repre¬ 
sents the results of lBrandt fc Huanel ( 2015 a) for the Prae- 
sepe open cluster. These authors interpret the spread as be¬ 
ing due to rotation, but show that if non-rotating isochrones 
are used, an age spread of ~ 200 Myr is in ferred. _ 

Fin ally, we show the results from iNiederhofer et ahl 
d2015al ) as upside-down (purple) triangles. While those au¬ 
thors did not interpret the derived age spreads as real, as 
they were largely interested in testing if such clusters had 
age spreads of 100s of Myr within them, they did find that 
the measured spreads were larger than expected due to pho¬ 
tometric errors. We took their results and de-convolved the 
total dispersion from the expected spread due to photomet¬ 
ric errors. Then we turned the resulting lcr dispersions to 
FWHM to more dire ctly compare with the other works. 
iBastian fc Silva-Villa (120131) also stud ied NGC 1856 with 
the same method as Niederhofer et ahl (l2015al ). namely in¬ 
vestigating the SFH using the MSTO and post-MS features. 
These authors found a dispersion of 35 Myr, correspond¬ 
ing to a FWHM of ~ 85 Myr. This is signi ficantly shorter 
than t hat found bv lMilone et al.l d 201 5l ) and ICorrenti et ahl 
(2015), with the apparent difference being due to the inclu¬ 
sion of the post-MS features, that show less variation than 
the MSTO. 

As can be seen, the observational data closely follow 
the prediction from the rotating isochrones. Hence, the data 
add quantitative confirmation for a key prediction of the 
rotating isochrones, that the inferred age spread is propor¬ 
tional to the age of the stellar system. This prediction can 
be further tested with high precision photometry of other 
young massive clusters in the LMC/SMC. For example, the 
slightly younger (~ 180 Myr) b ut si milar mass (~ 10 2 * * 5 M@) 
cluster, NGC 1866 dBastian fc Silva-Villal[2013l) should also 
show an eMSTO, but the extent of which should be less 
than NGC 1856 when translated into ’age space’. Another 
massive (> 10 5 Mg) young (~ 100 Myr) cluster that should 
show an eMSTO is NGC 1850, with the expected MSTO 
spread (when interpreted as an age spread) of 30 — 50 Myr. 

It is worthwhile noting that none of the clusters in this 
sample, or the lar ger sam ple of integrated spectroscopy ex¬ 
amined by IBastian et alJ (120131) , show any evidence of on¬ 
going star-formation. Hence, if the spreads were real, all 
star-formation must have stopped recently, with no clusters 
caught in the act of forming new stars. 


4 CONCLUSIONS 

In this work we used the Syclist stellar evolution mod¬ 
els to investigate the effect of stellar rotation on CMDs of 
star clusters at various ages. Comparing isochrones of non¬ 
rotating and rotating stars at different ages between 50 Myr 
and 1 Gyr we find that rotating stars, which have longer 
MS lifetimes due to enhanced internal mixing, can be mis¬ 
interpreted as younger stars. At the location of the MSTO, 

2 ICorrenti et al.l (120151) prefer a two-component single stel¬ 

lar population m odel, in order to be consistent with the 
iGoudfrooii et al] d2014l ) interpretation of the intermediate age 

clusters . This, however, is in conflict with the interpretation by 

Milo ne e t al. ll2015l) 


rotating isochrones closely follow younger non-rotating ones, 
causing a spread in the MSTO. This difference is propor¬ 
tional to the age of the isochrone in the considered age inter¬ 
val. In the light of the age spread interpretation, this would 
result in an increasing age spread with increasing age of a 
cluster. We compared literature measurements of apparent 
age spreads in young and intermediate age clusters with the 
predictions from the Syclist models. We found a strong 
correlation between the age of the clusters and the reported 
spread in ages, in close agreement with model expectations. 

We also investigated the effect that rotation has on 
the post-MS features of clusters and found that rotation 
produces evolved stellar tracks that are more compact in 
the CMD than the ones produced by the corresponding 
age spreads, although the MSTO region is similar in both 
cases, in agreement with the results i n ferred from young 
and intermediate age cl usters dBastian fc Silva- Villal 120131 : 
INiederhofer et aDl2015bl ). 

The model presented in this work that relates the eM¬ 
STO feature to rotation makes clear predictions that can 
observationally be tested. If rotation is causing the spread 
in the MSTO, also two young clusters NGC 1850 and 
NGC 1866 should show an extended MSTO, as well. How¬ 
ever, this extent should be smaller than the one found in 
NGC 1856, and of the order of 20 — 100 Myr, if interpreted 
as an age spread. Moreover, our model predicts different 
post-MS morphologies in the age spread and rotation in¬ 
terpretation, making the evolved stellar part of the CMD a 
valuable tool to discriminate between the two scenarios. 
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